General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



"Researches Into Planetary Atnospheres and Interiors" 


o 

o 

10 

1 “ 

(0 o> 

1 


A 

u o\ 


B "1 

as 

3 O 


Grant Number: NGL 22-007-2?.8 ^ 

o> 

m 

U 


SEMI-A14NUAL STATUS REPORT ?S 

Period: April 1, 1975 to September 30, 1975 


Submitted to: NASA Scientific and Technical Information Pacility 


P. 0. Box 8757 


Baltimore/V’ashinp.ton International Airport 
Baltimore, Maryland 21240 


Richard Goody, Mallinckrodfr' 



Professor of Planetary Physics 



w 

U4 


a 

ro 

03 


o 

« 

•a 


H 


4-1 

VM 

to 

u 

<« 

•o 

•-0 

■S 

•-( u 

-.J 

(-1 

<0 « 


tu 

9 > 



6 t-i 


O 

S « 


e- 

« = 



•H 


►-» 

B 



tp in 


to 

to 


M 

in 


as 



U 

to 


03 

04 • 


•« 

t_> U4 


lu 

K-f u; 


in 

04 to 

ca 

t»4 

M4 

Ul 

04 

H O 

• 


as m 

ro 


M 

to 


1 



(-) 

u 

IT) 

as • 

X 

CO 

•4 M 


U) 

04 

04 

•a 

Vi <4 


T“ 

04 

04 

1 

PS r- 

04 

04 

04 


u 

X * 


1 

04 *> 


■4 

CO M 

• 

CO 

O O 

> 

n 

X U. 

i*H 

as 

t-» O 

a 

4^ 

' 4« 03 

s 


Division of Engineering and Applied Physics 


Harvard University 

Cambridge, Massachusetts 02133 




I'Vi 






Cfej 

.fcSl 










October 1, 1975 



1. Introduction 


During the reporting period we worked towards a conclusion for two of our 
three principle tasks. !>r. Michael Flasar completed work on the transport of 
water through the martian boundary layer, and has moved to Goddard Space Flight 
Center. Alex Woronow completed his work on martian cratering and his Ph.D., 
and has now moved to the Lunar and Planetary Laboratory of the University of 
Arizona. 

The principle Investigator also completed an Investigation of the photometric 
properties of lo in collaboration with Dr. Yuk Yung. 

The third major branch of our work, planetary observing, will greatly increase 
In scope in future j-ars. We also intend to Investigate further the problem of 
water transport on !!ars by collaborating with the Viking project. 

2-1 . Echelle spectrograph and Kron Camera 
(Dr. David Latham) 

Good progress was made over the past six months on the project to develop 
instrumentation for Planetary Spectroscopy using the 61-lnch telescope at 
Agassiz Station. Tlie echelle spectrograph vas assembled over the summer and put 
into operation on the telescope in October. A few dozen good spectra were 

O O 

taken, including some plates of lo covering 4000 to 7000 A at about 2 A/raci, 

The Kron Camera \;as also assembled over the simmer, and an S~ll cathode with 
about 13% quantum efficiency v/as applied successfully in August. Use of the 
Kron Camera at the telescope is waiting for completion of the 30KV voltage 
divider and delivery of some special high voltage cables. 

The source of the distortion (about 12%) in the Mt. Hopkins Kron Camera was 
traced to the corrector lens, which was originally designed at the Naval Ob- 
servatory to map the focal plane of their 40-inch Ritchey-Chretien telescope 
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onto the curved cathode. This required a t»\' 0 -elenent lens with a wild aspheric 
in order to correct for the astlf^atisra of the 4C-inch. A new single-element 
lens has been designed specifically to match the camera to the spectrograph. 

The lens has been fabricated and will be tested shortly. This represents a 
major breakthrough if it works out, because it v;ill very much simplify the 
tracing and reduction of echelle spectra taken with the Kron Camera. 

2-2 . Properties of the martian surface 
(Alexander Woronow) 

Work has been completed on two important aspects of the analysis of martian 
crater statistics. First, craters in the diameter range 10-250 km were demon- 
strated to have a log-normal size-frequency distribution, not a log- log dis- 
tribution as previously assumed. This implies that dates derived from straight 
isochrons on a log-log model using craters of less than 40 km diameter could 
be systematically underestimated by an order of magnitude. Second, using the 
log-normal model, regions on Mars were sought that had statitstically anomalous 
crater size-frequency distributions. Tixelve such regions were located. Most 
important among these are 1) the northeim volcaiic plains; 2) the area south 
of -30“ latitude; and 3) a region centered on +35“ latitude, 10“ west longitude. 

The first of these, the northern volcanic plains, has an over- abundance of 
small craters. Tills is most likely caused by a decrease in the mean radius of 
the impacting bodies with time. This is the first evidence for such an effect. 

The second region, that area south of -30“ latitude, has a deficiency of 
small craters. This result is consistent with the observations of Soderblon, 
et al. , that an eollan blanket covered many craters up to 10 km diameter. This 
study, hov7evcr, demonstrates that some regions have had craters beyond 20 km 
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diameter affected. 

The third region, that centered on +35* latitude, 10® west longitude, dis- 
plays a deficiency of large craters. This observation, coupled with associated 
morphologies of channels and knobby terrain, makes this region the most likely 
site fotjnd for large quantities of subsurface ice. 

These and additional results were incorporated in the Ph.D. thesis of Alex- 
ander Woronow, which was defended this September. A condensation of the thesis 
is in preparation for submittal to JGR . Alex Woronow has accepted a post-doctoral 
position with Robert Strom at the Lunar and Planetary Laboratory of the University 
of Arizona, where he is currently examining the implications of crater saturation 
and equilibrium through computer gener;'.teJ stochasitc model, and examining the 
time evolution of icipacting bodies as recorded on Mercury. 

Ho further work on this part of the project will take place at Harvard. 

2-3. Martian water 

(Dr. Michael Flasar, Ralph Kahn) 

A paper on the diurnal behavior of Martian water, described in the previous 
report, has been accepted for publication by Planetary and Space Science . 

As of October, Dr. Flasar has taken a position at the Goddard Space Flight 
Center v?ith R. Hanel's g' oup as a national Reserach Council Resident Research 
Associate. He is studying the problem of thermal convection within the interior 
of Jupiter. Special emphasis is being given to the question of how the planet's 
rapid rotation affects the convection flow; in particular, does rotation induce 
a meridional variation in Jupiter's flux of internal heat? The effects of pos- 
sible first-order phase transitions in inhibiting thermal convection v;lll also 
be investigated. Dr. Flasar' s work was begun in August and September under 
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this grant and any relevant publications ^/lll acknowledge this support. 

Further work along two main lines is proposed in this area. Palph Kahn, 
as part of his thesis research, will Investigate climatological tine scale 
transport of water, basrd upon the understanding of the boundary layer reached 
by Flasar. 

As an Irrr’ediate problem, Kahn and Goody have joined the Vih.lng imaging tet a 
as guest investigators. This results from the prediction of our previous 
work that ice mists should exist between about 2400k and IHOOk, even in tropical 
regions. Enquiries were made through Or. Soffen and Professor ''uteb as to whether 
the V'iking camera was suitably programmed to tecerd these mists — if they exist. 

We have nov; entered into a collaboration with Er. James Pollack of Ames P.e- 
search Center. We will jointly design observing sequences and the programs 
needed to invert the data. Ralph Kahn will spend most of the summer at JPL working 
with the Viking team. 

Ralph Kahn’s salary is being paid from fellowship funds. Tlic Viking project 
has, hW'/ever, no money for his expenses, and these \7ill be charged to KGL 22- 
007-22S. A visit by Goody and Kahn to Ames P.esearch Center has already been 
made with support from the grant. 

2-4. Photometric properties of the surface To and their influence on line 
formation in the atmosphere 
(Dr. Yuk Vung) 

The following is the abstract of a paper x^rltten by Yuk L. Yung and Richard 
M. Goody, and submitted for publication in Icarus . 

Abstract 


Wc give a quantitative theory of line formation in an atmosphere above a 
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surface with backscattcring properties. Sufficiently high spatial and spectral 
resolution spectra of resonance lines in To region A can yield data on the sur- 
face scattering properties as well as the nunber density of scattering molecules. 

We discuss nacroscoplcally homogeneous models of scattering from the sur- 
face of lo and conclude that multiple reflection from crystal facets is the most 
likely cause for the observed phase variations of the geometric albedo. 

This work is important for the interpretation of spectrt as discussed in 
section 3, below. 

2-5. Imaging of lo in the NaD lines 

(Frank Hurcray) 

Our picture of the sodium system is now similar to the hydrogen torus 
proposed by McDonough and Brice for Titan. Its form will depend upon the 
escape process from lo and the forces encountered subsequent to escape. 

We cannot predict the form or behaviour of the forces “»ith any confidence, 
and the one-dimens ior.al data which have been obtained by Wehinger and 
Wyckoff (1974' ■'eckler et al . (1974, 1975) and by us are tantalizingly 
inadequate for visualizing the spatial distribution of sodium. Two- 
dimensional images from which we can study the time dependence and 
morphology of the sodiim system within 2 to 3 arc -minutes from lo are 

essential. Unfortunately they are very difficult to obtain. 

0 . 

We have taken images through a 7A filter centered on the D 2 line, 
using the Mt. Hopkins Kron camera. We believe that we cam just detect 
the sodium emission near to lo but light scattered in the telescope from 
Jupiter is very strong even when we take precautions to eliminate it. 

The Na line can be picked out with a spectrograph because the high 
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resolution discriminates against the continuum radiation (the lo lines are 
less than 100mA wide) but, more importantly, the line and adjacent 
continuum are recorded side by side and simultaneously on a very small 
area of the detector. We could probably get a good sodixim image with a 7 a 
filter if a continuum image could be interwoven with it and differences 
for ted at a later stage. This is a general problem in the imaging of 
dif:fuse astronomical objects. 

We have designed and are building a two-dimensional imaging system 
which has many of the advantages of a high-resolution cpectrometer. 

Consider a conventional spectrometer whose entrance slit is replaced 
by a grid of clear and opaque strips perpendicular to the dispersion, and 
spaced by a distance, A . If this grid were ill’amiuated with strictly 
monochromatic light and a ceunera placed at the exit aperture the grid would 
be imaged exactly and a one -dimensional scan along the dispersion would 
produce a record such as is shown in Figure 1(c), This is the case e -*• “ 
where e is the ratio of A to the width of the line (see Figure 2.). 

Now consider a wide line ( C small) with the same central intensity. 

A scan of the image will look like Figure 1 (b) . The average signal must 
increase proportionally to the line width of the illuminating line but, for 
1 » ethe ripple on the record will rapidly decrease because of overlap. 
Provided that we have a linear deteccor with large dynamic range (as is the 
case for the Kron camera) we can concentrate attention on the ripple. 

For example, it would be possible to convolute the record with a sampling 
function such as that shown in Figure 1 (c) , which has a zero average, and 
therefore no response to a continuum. Figure 2 shows the result of scanning 
such a function over the record in a computer simulation. 
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Bxe important feature of Figure 2 is the rapid fall off near e ■ 1. 

If we pick A so that e « 4 for the Na lines there is a dramatic discrimination 

-2 

against wider lines. 1116 secondary peak near c «■ S x 10 should not be 
realized if an appropriate blocking filter is used (about lOA wide will 
suffice if we optimize for O.lX lines) . The slow fall off of response as 
e increases is not important since we do not anticipate lines narrower than 
the NaO lines. The problem is to discriminate lo emission against solar 
lines which arc at least 10 times wider. 

Recovery of the image may follow a process similar to that described 
if performed on the computer. An alternative # optical procedure is 
illustrated in Figure 3. 

We are currently building the equipment needed for this imaging device. 

Dr. David Latham of SAO is constructing an improved version of the 
Mt, Hopkins echelle spectrograph for the Agassiz telescope. In collaboration 
with R. Brown he is also constructing a Kron camera for use on the 
spectrograph and for other piurposes. It is hoped that both instruments wi.' 
be ready by the end of the svimmer. 

While they are being constructed we have designed and constructed 
entrance grid to replace the entrance slit system. We have also built an 
offset guidar. This work is of interest to members of the SAO because of 
its application to the general problem of nebula imaging. Dr. Lawrence Mertz 
has constructed an alternative system based upon a wedge interferometer. 

We expect to work both projects on a collaborative basis. 

Since we are not dealing with very faint or small objects the sky 
conditions in Massachusetts should not be as disadvantageous for this work 
as it is for more conventional astronomy. We hope that much of our observation 
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can therefore take place at the Agassiz station although we will probably 
need to take advantage of the better skies of Mt I’opkins at a later stage. 

In recent months our v;ork has advanced rapidly. Dr. Mertz's interferometer 
has been tested by Mertz and Murcray on the Agassiz 61" telescope. Murcrav has 
taken the instrument to Mt. Hopkins and has used it with an inage tube. Fringes 
have been obtained of the Crion nebula and also in the general area of lo. Plates 
are nov; being processed. The laser analysis device has been constructed in the 
laboratory. 

V7e are still some way from having solved all of the probler..- associated v?ith 
the Interferometric inage. Hov7cver, v?e are solving then one by one and perfor- 
mance seems to correspond to prediction at each stage. 

Simultaneously with our work on the interferometer, the imaging system de- 
signed for the echelle spectrograph is nearing completion. The offset guide is 
complete and the entrance grid has been purchased. A three-stage image tube 
belonging to SAG would be ideal as a detector. Fork has to be done on the 
image tube, hov»ever, before it is ready for use. ’’c are also still vraiting for 
final tests on the Kron tube. 

It is difficult to predict where this project ^.ill be in si?: months or \;hich 
of the two instruments and alternative detectors will be m.ore suitable. V'e are 
advancing on several fronts at the same tire and the chances of success look good. 

3. Future vrork on planetary observing 

With the agreement of hr. Robert Fellows, the work, of Dr. Robert Drown “fs 
being transferred to KGL 22-007-228. The three main objectivos of this w’o t,re 


described belov/. 
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3- 1 , Search for liuorcscent cn'lsslons 

The sodium emission from lo and the Lyman-alpha emission subsequently 
discovered by Caurlson and Judge (1974) during the Pioneer 10 flyby both 
draw attention to the possibility that other bodies ir the outer sol 2 u: 
system may have tenuous atmospheres. If Jupiter satellites are icy, as 
prc'p>sed by Lewis, the ice is probably evaporating at the present time 
and tie satellites should oo surrounded by escaping atoms which could 
fluoxosce. The same may also be txnae of seme of the asteroids. 

We propose to search for fluorescent emissions from icy objects, 
satellites or asteroids. The species Li, Ca, K, Mg, 0 have already been 
sought on lo, with negative results. We also have some spectra of the other 
three Gallilean satellites which we are in process of comparing. 

The equipment available to us is ideal for exploratory spectroscopy. 

A larger telescope than 60" would be valuable for very faint objects but 
there are many for which we can obtain good spectra with a reasonable 
exposure time. Cassegrain echelle spectrographs are fast and more effective 
than spectrographs on many larger instruments. The combination of 
cross-dispersion and a Kron camera allows us to store up to 1000 A of 
spectrum with photometric accuracy and a dynamic range of lO^. 

The problem with this system is the immense amotont of high quality 
data stored on a single Kron plate. It is too good to be used only for 
visual inspection and we have therefore developed techniques for conputer 
processing. 

We have modified the Mann microdensitometer of the SAG so that spectral 
scans can be put automatically onto tape. Programs have been written 
which can translate and scale ^e data so that two spectra from the same 


10 


re9ion can be registered to an accviracy higher than the spectral resolution. 
Two spectra from different objects cam then be differenced. A section of 
two spectra and the difference spectrum arc shown in Figure 4. 

With these automatic data handling techniques we will be able to 
analyze a large number of spectra with minimum effort on our part, and 
we propose, over the next two years to e:diaust the spectral range of the 
present Kron tube cathode for Jupiter satellites and the larger asteroids, 
looking both at the discs and the nearby sky. 

V2 . Ro tation rates of T»ranus ar.H 

The characterization of dynamical regimes in a planetary circulacion 
depends upon the rotation rate (Stone, 1973) . In the case of Uranus and 
Neptune surface markings cannot be seen sufficiently well to determine the 
rotation rate and there is considerable doubt as to the significance of tlie 
published data. 

The problem is reviewed by Newburn.. and Gulkis (1973). Data from time 
variability of the light curve of Uranus are contradictory. The Doppler shift 
work of Moore and Kenzel was performed in 1928 and 191C v;ith prism 
spectrographs. It is very hard to believe that tlieir measurements are of 
any significance although they claim 5 to 7% accuracy. 

We could make very much more accurate measurements with modern 
spectrographs. With well resolved discs 1% accuracy should present no 
difficulty. Uranus can just be resolved under normal conditions, but the 
disc of Neptune cannot be explored. 

The least possible distortion by' the spectrograph along the slit 
is required to take advantage of disc resolution and the best instrument 
known to us is the stigraatic echelle spectrograph of the KPNO 150". We 
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•hall seek time on this instrument in Spring, 1976. Preliminary measurements 
will be made at Mt. Hopkins in May of this year. 

An edtemative approach, which v;e shall probably have to use for 
Neptune, is to combine all the .’’^ghi; from the planet and to fit the 
resultant line shape nmerically. To do this we must know the direction 
of the rotation axis. It is commonly assumed that the rotation axis is 
■ perpendicular to the plane of the satellites. Unfcrtunatelv. the only data 
capable of confirming this, the direction of faint bands on Uranus, fails 
to do so. Unless observations can be made from the proposed LST many 
accurate spectrometric observations over several decades will be needed to 
solve this problem. It is time to start the process. 

3-?. Thermal anornljes or» Vonttc; 

The dense atmosphere, the ubiquitous clouds and the slow rotation rate 
of the planet put the Venus atmosjiiere in an interestingly different para- 
metric range from Earth (Gierasch, Goody and Stone, 1970) . This has led to 
a series of studies of possible Hadley circulations in the lower atmosphere 
(Goody and Robinson,, 1966; Stone, 1974; de Pivas, 1975), The proposed 
Pioneer-Venus missions of 1978 have been designed to explore this deep 
circulation. There is little further than can be done to investigate this 
phenomeipn f rom Earth before this mission is completed, 

A those levels where solar radiation is absorbed (60 km above the 
surface) a diffuse circulation, having a boundary-layer character, can be 
observed. This is the remarkable 4-day circulation (Young and Schubert, 
1973), about which there are many theories: eventually, of course, a complete 
theory must embrace both the boundary layer and the interior, but such 
%#ork is still in its infancy. 
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Any doubts as to ths reality of the 4-day circulation were put to rest 
by the Mariner 10 fly-by (Murray ejt aju 1974) . Although much interesting 
detail was revealed in the blue images they do not give us any obvious hints 
as to the physical causes oi this circulation. 

The 4-day circulation was first detected from Earth, although the 

observations were difficult to interpret. This emphasizes, howeveo, the 

need to pursue all available observatlcnal data, particularly v^en the 

phenomena appear to be time-dependent emd therefore related to dynamical 

effects. One class of data which shows this promise, but has not been 
fully exploited, is represented by the thermed maps of Murray, Wildey and 

Westphal (1963) and Westphal, Wildey and Murray (1965). 

Attempts to analyze these maps in some detail have been made by Goody 
(1965) and Ing -rsoll and Orton (1974) . Limb darkening, solar related 
effects and equator- to -pole gradients are fairly well described observationally 
although theoretical explanations are lacking. 

Features of these maps which have been little investigated are the 
thermal anomalies, such as that seen near the S-pole by Murray et a^. (1963). 

In view of the paucity of data on Venus and the probadsle relevance of any 
large scale temporal variations to the planetary circulation, either as 
cause or effect, we believe that these phenomena should be further 
investigated. It has not been clear, for exan 5 >lc, whether thermal anomalies 
defirltoly exist. If they do, we may be able to understand their relationship 
to the circulation if we know their size, their intensity, their lifetimes 
and their preferred motions, if any. As pointed out by Ingersoll and Orton 
we cannot distinguish, from measurenents at a single wavelength, whether 
the effect is due to c -.ges in cloud density or changes in cloud temperature. 
Measurements at more than one wavelength would be required for this purpose. 
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fii >d as a conseqaence we have concentrated attention on "black” detectors which 
canf in principle, be used at any wavelength. 

Apart from the decision with respect to the detector o\ir other two 
major decisions have been to use the Mt. Hopkins 60* telescope because of its 
availability and to attenpt image stabilization- for accurate and reproducible 
mapidng. The large diffraction image {approximately 2" diameter) of the Mt. 
Hopkins telescope is unsatisfactory, but we judged ready availability to be 
more inportant than hi^ resolution when studying time-varying pheromena. 

Ihe decision to stabilize the image led to the use of the coude focus. 
He had hoped that the large signal from Venus would allow us to compensate 
for other difficulties associated with this focus. For technical reasons 
it led us into a slow scanning mode with an undesiredsle amount of the 
1/f detector noise which gives image "streakiness”, despite careful image 
processing. This hats set the limit of accuracy to our measurenents. 

Our equipaent consists of a Low-type bolometer with an image stabij-izer 
and computer-controlled raster scanning and chopping equipment. A schematic 
diagram is shown in Figure 5. We obtained our first series of thermal 
maps in the period Novej?i>er 1973 to January 1974. Since that time we have 
been developing the many software programs needed to remove instrumental 
distortion, reduce the noise, remove liznb darkening tto a first approximation) 
and to laap onto an oriented roercator projection. 

Figure 6 shows one example of a thermal map. Note that most of the 
limb darkening has been removed amd the remaining energy contours are close. 
Well-known, solcir-related effects can be found in this image, but there 
is also a general warming in the S mid-latitudes which is reproducible .*nd 
has not previously been reported. The anomaly in the north-west segment also 
repeats in neighboring records but only for a few days. We have therefore 
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confirmed the existence of tine-dependent anomalies and have sho%m that 
the 60" telescope can resolve at least the larger ones. 

During the next grant period we have both short -tern and long-term 
objectives. Our existing equipment r being based on conputer instructions, 

C£m be used in many soodes. This summer a scanning secondary will be added 
to the Mt. Hopkins telescope which can be controlled with our equipment 
to perform raster scans. The detector can then be used at Cassegrain. Image 
6 tabilization will not be possible but hopefully scanning can be rapid 
enough to make this unnecessary. These procedures should give great improvements 
in signal-to-noise ratio and the rate of acquisition of data. 

On the longer term we would like to use the Mt. Hopkins Multi-Mirror 
Telescope, provided that the images can be phased in the infrared. If so, 
resolution as good as 1 arc-sec could be available for thermal imaging. 
Approximately one year from now we plan to become involved in the development 
of the infrared imaging system for this telescope. We do not intend to 
contribute hardware, however, since budgets are already set up for this purpose. 
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Figure 4. lo and Europa spectra in the same wavelength range and the 
difference spectrum after processing. 


Figure 5. Thermal mapping equipment. 
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